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Lidar  detected  spike  returns 


William  P.  Hooper  and  Glendon  M.  Frick 

Naval  Research  Laboratory,  Remote  Sensing  Division,  Washington  DC  20735,  USA 

bill.hooper@nrl.navy.mil 

Abstract.  Lidar  measurements  of  spike  returns  from  clear  air  are  presented.  These  spikes 
occur  infrequently  (approximately  one  in  hundred  returns)  but  provide  returns  that  are 
significantly  stronger  (occasionally  an  order  of  magnitude  larger)  than  the  average  aerosol 
backscatter  signal.  The  spike  density  is  5.7e-3  spikes  m‘3  for  backscattering  cross  sections 
estimated  to  be  between  0.003  and  0.080  mm2  sr'1.  A  modified  form  of  the  lidar  equation 
which  includes  returns  from  large  particulates  is  presented  and  the  probability  distribution  for 
the  spike  magnitudes  is  derived  from  five  million  measurements. 

Keywords:  lidar,  near-infrared,  insect,  large  aerosols,  ultra-large  aerosols,  spike  returns. 


1  INTRODUCTION 

Strong  isolated  returns  have  been  noted  from  the  earliest  days  of  active  remote  sensing.  In  the 
1940's  when  radars  were  being  developed,  strong  "spike"  returns  from  clear  air,  where  no 
return  was  expected,  were  whimsically  called  "angels"  [1].  Measurements  showed  that  the 
characteristics  of  the  angels  depended  on  the  radar  transmit  wavelength  and  were  caused  by 
hard  target  returns  or  "spot  angels"  from  birds  [2],  insects  [3],  and  more  diffuse  returns  from 
refractive  structures  [1]. 

Early  lidars  also  occasionally  detected  spikes  from  birds,  insects,  hydrometeors  and  ice 
particles;  however,  since  much  of  the  research  focused  on  volumetric  scattering  from  large 
numbers  of  smaller  aerosols  which  impact  atmospheric  optics,  these  strong  returns  were 
normally  not  reported.  Unlike  radar,  lidar  does  not  detect  scattering  from  refractive  structures 
[4]  and  thus  the  lidar  returns  are  caused  by  either  molecules  or  aerosol  particles.  (In  this 
paper,  the  term  "aerosols"  will  denote  small  atmospheric  particles  and  the  term  "particulates" 
will  denote  larger  particles,  insects  and  birds). 

The  spike  returns  discussed  in  this  paper  are  strong  isolated  returns  detected  by  a  short- 
pulse,  incoherent  lidar.  Figure  1  shows  examples  of  the  largest  and  smallest  spikes.  These 
isolated  spike  returns  were  detected  in  clear  air  while  testing  an  eye-safe  lidar.  The  spike 
returns  are  significantly  stronger  than  background  scattering  and  occurred  with  no  visually 
apparent  scattering  source.  Initially  concerns  were  that  the  spike  returns  could  be  caused  by 
instrumentation  and  the  detector,  detector  power  supplies,  and  the  digitzer  were  tested  and 
yielded  no  instrument  problems  that  could  generate  spike  returns.  Since  the  spikes  occurred  at 
the  same  ranges  that  returns  are  seen  from  atmospheric  aerosol  structures,  the  larger  spikes 
are  apparently  caused  by  scattering  from  one  or  more  "sub-visible"  particulate  targets  and 
contrast  with  the  smaller  backscatter  signals  from  nearby  measurement  volumes  that  only 
have  the  smaller  aerosols  and  have  no  large  particulates. 

In  this  paper,  we  will  discuss  the  strength  of  these  returns  and  their  probability  of 
occurrence.  Section  2  discusses  the  lidar  including  features  which  may  explain  the  spike 
detection.  Section  3  discusses  a  modified  lidar  equation  and  the  data  analysis  technique  used 
to  characterize  the  spike  returns.  Section  4  discusses  the  measurements  and  finally  section  5 
summarizes  the  measurements  results. 
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Fig.  1.  Lidar  profiles  showing  the  largest  and  smallest  spike  returns.  The  background  signal  in  the 
blue  curve  is  similar  to  the  background  signal  in  the  red  curve;  however,  the  blue  curve  was  scaled 
to  offset  it  from  the  red  curve  and  allow  both  curves  to  be  plot  together. 


System  Block  Diagram 


Fig.  2.  Schematic  of  the  NRL  Eye-safe  Volume  Imaging  Lidar  (NEVIL)  System. 

2  INSTRUMENT  DESCRIPTION 

An  existing  scanning  lidar,  which  has  been  used  in  a  series  of  field  experiments  (for  example, 
see  Hooper  and  Martin,  1999)  [5],  has  been  modified  to  form  an  eye-safe  lidar.  The  original, 
non-eye-safe  transmitter  used  a  Nd:YAG  laser  to  transmit  at  1064  nm.  The  new  eye-safe 
system  is  build  around  coupled  deuterium  cells  [6]  that  shift  the  optical  pulses  at  1064  nm  to 
1560  nm.  Figure  2  shows  the  schematic  and  Table  1  lists  the  parameters  for  the  NRL  Eye-safe 
Volume  Imaging  Lidar  (NEVIL)  system.  The  Nd:YAG  laser  pulse  is  split  between  two  cells: 
thirty  percent  is  directed  into  the  seed  cell  with  the  remaining  energy  directed  into  the  pump 
cell.  A  backward-propagating  beam  is  generated  by  the  seed  cell,  amplified  and  combined 
with  the  pulse  from  the  pump  cell  which  provides  a  light  pulse  that  should  replicate  the 
divergence  of  the  original  Nd:YAG  laser  beam.  The  Nd:YAG  laser  is  injection  seeded,  which 
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reduces  the  laser  line  width;  a  narrow  line  width  is  essential  and  allows  the  Raman  conversion 
processes  to  have  an  operating  efficiency  of  20%  or  better.  When  we  measure  the  beam 
quality,  we  find  a  beam  divergence  of  100  (trad  or  less,  which  differs  from  the  value  that  we 
expect.  It  is  possible  that  the  Raman  conversion  process  is  more  efficient  for  the  more  highly 
collimated  photons  and  thus  changes  the  divergence  of  the  output  beam.  Since  the  pulses 
counter  propagate,  the  output  pulse  is  shorter  (3  ns)  than  the  original  Nd:YAG  pulse  (8  ns).  A 
Galilean  beam  expander  increases  the  beam  diameter  from  approximately  2  cm  to  6  cm.  The 
resulting  1560  nm  beam  is  transmitted  into  the  atmosphere. 


Table  L  NRL  Eye-safe  Volume  Imaging  Lidar  (NEVIL)  System  Parameters 


Transmitter 

Type 

Raman-shifted  Nd:YAG 

Pump  laser 

Continuum  8010  (injection  seeded) 

Pump  linewidth 

0.003  cm1 

Output  Wavelength 

1560  nm 

Pulse  length 

3  ns 

Energy 

200  mJ/pulse 

Repetition  Rate 

10  Hz 

Divergence  (for  output  beam) 

100  |irad  (full  width  at  half  maximum) 

Receiver 

Telescope 

Cassegrainian  (diameter  0.36  m,  focal  ratio  f/11) 

Detector 

InGaAs  APD,  diameter  200  pm 

Detector  Collecting  Optics 

Drum  Lens  (BK-7,  4  mm  diameter,  curvature  leading 
surface  2.5  mm) 

Field  of  View 

100  pirad 

Digitizer 

50  MHz,  12  bit 

The  light  backscattered  from  the  atmosphere  reflects  off  a  dual  mirror  scanning  system,  is 
collected  by  a  0.3  m  Cassegrainian  telescope  and  the  backscattered  signal  is  focused  onto  a 
200  pm  diameter  InGaAs  avalanche  photodiode  (APD).  A  small  drum  lens  is  mounted 
directly  above  the  APD  and  used  to  increase  the  APD  collection  efficiency.  (Drum  lens  are 
machined  from  ball  lens  into  a  cylinders  and  are  easier  to  align  with  cylindrical  detector 
packages.)  The  receiver  field  of  view  is  approximately  100  (trad. 

The  lidar  system  was  tested  on  a  400  m  range.  During  these  tests,  we  found  that  the  beam 
divergence  was  also  approximately  100  (trad  and  that  the  energy  distribution  across  the  beam 
has  a  Gaussian  distribution.  The  transmit  beam  fully  illuminates  the  receiver  field  of  view; 
thus,  while  the  backscatter  signal  from  particulates  from  the  center  of  measurement  volume  is 
stronger  than  the  signal  from  the  side  edges  of  the  measurement  volume,  both  signals  should 
be  detectable  if  the  particulate  backscatter  is  significantly  stronger  than  the  background 
aerosol  signal.  Therefore,  we  assume  that  the  measurement  volume  (nr2  (Acp)2  Ar)  is 
determined  by  the  range  (r),  the  receive  the  field  of  view,  where  the  angle  (Acp- 50  |_irad)  is 
half  the  receiver  field  of  view,  and  the  range  resolution  (  Ar). 

Prior  to  our  measurement  runs,  a  green  laser  pointer  is  placed  in  the  far  field  and  is  aimed  at 
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the  lidar.  The  green  laser  is  focused  both  onto  the  InGaAS  APD  detector  and  into  the  D2  cells 
of  the  transmitter.  With  the  pump  cell  removed,  a  pinhole  is  placed  at  the  location  of  the  focus 
spot  for  the  1064  nm  YAG  beam.  When  the  focused  green  beam  is  focused  both  on  the 
detector  and  the  center  of  the  Raman  cell,  the  system  is  aligned.  This  alignment  maximizes 
the  overlap  function  of  the  transmitter  and  receiver. 

The  new,  eye-safe  lidar  differs  from  the  old  non-eye-safe  system  in  that:  the  wavelength  is 
shifted  from  1064  nm  to  1560  nm,  the  pulse  length  is  reduced  from  8  ns  to  3  ns,  beam 
divergence  is  reduced  from  1  milliradian  to  100  |irad,  and  range  resolution  is  reduced  from 
7.5  m  to  1.5  m.  These  changes  allow  the  detection  of  small  aerosol  structures,  such  as  those 
generated  by  breaking  waves.  However  as  we  explain  in  the  following  two  sections,  these 
changes  also  make  the  spike  returns  from  large  particulates  significantly  more  detectable. 

3  THEORY 


Lidars  operating  in  the  near-infrared  receive  backscatter  signals  from  molecules  and  aerosols. 
All  the  measurements  presented  in  this  paper  made  within  300  m  of  the  surface  and  are  within 
the  Earth's  boundary  layer  where  the  aerosol  signal  is  especially  strong.  The  lidar  equation  for 
volumetric  aerosol  scattering  is  [7]: 


P(r)=Kr  2 PT[r)Arexp 


J  E(r')dr' 

r'=0 


+n  r 


(1) 


where  P  is  power  return  for  volume  backscatter,  K  contains  system  constants  and  laser  output 
energy,  pT  is  total  volume  backscatter  coefficient  from  all  aerosols  and  particulates,  e  is  the 
volume  extinction  coefficient,  and  n  is  the  offset  caused  by  sky  and  electronic  noise.  The 
volume  backscatter  and  extinction  coefficient  are  impacted  by  both  aerosol  and  molecules; 
however  in  the  atmospheric  boundary  layer,  the  molecular  backscatter  and  extinction  in  the 
near-infrared  are  typically  more  than  a  hundred  times  weaker  than  the  aerosol  backscatter  and 
extinction;  thus,  the  molecular  scattering  can  be  neglected. 

The  density  of  aerosols  depends  strongly  on  the  size  of  the  particles.  In  the  atmosphere's 
clear  air  boundary  layer,  there  are  often  billions  of  small  particles  (less  than  fifty  microns  in 
size)  in  a  cubic  meter  and  a  very  small  number  of  particles  larger  than  fifty  microns  in  a  cubic 
meter  (in  fact  this  density  may  be  substantially  less  than  one).  The  volume  backscatter  can  be 
arbitrarily  split  between  the  volumetric  scattering  from  large  number  of  small  particles  and 
the  sum  of  scatters  from  a  small  number  the  large  particulates: 

N 

(Sr(r)=iS0(r)  +  n1r  2(A(p)~2(Ar)~1'^l  crni  ,  (2) 

1=1 


where  po  is  the  volume  backscatter  coefficient  from  smaller  particles  (aerosols),  i  is  the  index 
for  the  particulates,  N  is  the  total  number  of  particulates  in  the  measurement  volume,  and  oni 
is  the  backscattering  cross  section  for  ith  particulate  (with  units  of  mm2  sr1).  For  this 
calculation,  the  measurement  volume  is  assumed  to  be  uniformly  illuminated,  see  the 
discussion  in  Section  2.  The  backscattering  cross  section  is  defined  as  the  differential 
scattering  cross  section  [8]: 


a 


T T 


dcrT(n) 

dQ 


(3) 


where  oT  is  total  scattering  and  the  derivative  is  taken  for  backscattering  into  the  solid  angle 
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element  dQ. 

Using  Eq.  (2),  an  alternative  form  of  the  lidar  equation  is: 


'  N  \ 

1  r  \ 

P{r)=Kr  2 

PoAr+  ti1  r  ^lAcp)  2^  ani 
\  i=i  ) 

|eXp| 

—  2  J  elr'Jdr' 

i  r~°  / 

Equation  (4)  explains,  in  part,  why  the  new  eye-safe  lidar  detects  spike  returns  when  the  older 
system  did  not.  As  the  range  resolution  and  the  divergence  angle  decrease,  the  volume 
backscatter  term  (poAr)  decreases  and  the  particulate  scattering  term  (  n'V2(A(p)‘2 )  increases 
and  the  ratio  of  the  particulate  scattering  term  to  backscatter  term  is  about  five  hundred  times 
larger  for  the  new  eye-safe  lidar  than  for  the  older  lidar. 

If  there  is  only  a  single  large  particulate  in  a  measurement  volume,  the  backscattering  cross 
section  becomes: 


nr2  A(p2  P^Ar 


‘Ps[r)-PB[rA 


\ 


PB(r)-n 


where  Ps  denotes  the  power  return  from  both  a  spike  and  the  background  scattering,  PB 
denotes  the  power  return  the  background  scattering  only  and  n  is  noise  offset  derived  from 
average  return  values  beyond  5  km.  The  background  return  is  calculated  from: 


4(p 


+p  i 


*Bnr4r  n-ej \r  ]-s)+p[r]-4i  n-3/j  ; 

where  j  is  the  range  index  for  the  spike  return.  The  response  time  of  the  lidar  receiver  is 
limited  and  large  spike  returns  are  followed  by  a  short  decay  in  the  signal,  which  impacts  the 
next  range  measurement;  therefore,  to  avoid  biasing  the  data,  the  background  signal  is  derived 
from  data  which  occurs  before  the  spike.  The  magnitude  of  return  bins  after  spike  returns  is 
also  checked  to  ensure  that  the  return  decreases  to  values  which  are  close  to  the  pre-spike 
background.  In  addition,  the  resulting  spike  observations  are  filtered  to  remove  duplicate 
measurements  of  large  spike  returns  at  consecutive  range  bins  and  from  profiles  at  the 
subsequent  measurement  times  with  the  same  ranges. 

This  paper  uses  Eq.  5  to  estimate  the  backscattering  cross  section.  Most  of  the  variables  are 
determined  by  measurement  and  the  errors  associated  with  the  variables  are  known  and 
limited.  The  largest  error  is  associated  with  determining  a  value  for  the  volume  backscatter 
coefficient.  This  value  is  assumed  to  be  constant  and  is  selected  from  the  report  written  by 
Wright  et  al.  [9]  which  give  the  coefficient  as  a  function  of  wavelength  and  environmental 
conditions.  The  error  in  the  value  of  backscatter  coefficient  could  be  as  large  as  an  order  of 
magnitude.  As  a  result,  the  error  in  value  of  the  derived  backscatter  cross  section  could  also 
be  as  large  of  as  an  order  of  magnitude.  Variations  of  volume  backscattering  depend  strongly 
on  the  variations  relative  humidity  [10].  From  the  estimated  humidity  variations  during  the 
measurement  period,  the  errors  caused  by  variations  in  volume  backscatter  coefficient  are 
expected  to  be  less  that  a  factor  of  two. 


+P|r 


4  MEASUREMENT 


On  17  April  2008,  lidar  observations,  which  were  by  NRL  Eye-safe  Volume  Imaging  Lidar 
system  made  over  the  Chesapeake  Bay  in  a  clear  air  boundary  layer,  detected  1145  spikes  in 
20  minutes  at  ranges  between  1.5  and  2.25  km.  The  winds  on  this  day  reported  by  NOAA  at 
Thomas  Light  (which  is  twenty  kilometers  North  of  lidar)  were  approximately  4  m/s  from  the 
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South  and  the  wind  direction  was  almost  perpendicular  to  the  trajectory  of  the  lidar  transmit 
beam.  Since  the  lidar  operate  at  10  Hz,  the  wind  can  move  the  aerosols  from  one  lidar 
measurement  volume  before  the  next  lidar  measurement  is  made.  During  these  measurements 
the  lidar  was  aimed  at  fixed  elevation  and  azimuth.  Figure  1  shows  spikes  typical  for  this  day, 
which  are  also  typical  of  the  "clear-air"  spikes  detected  by  the  eye-safe  lidar  over  the  last  few 
years.  While  we  have  observed  hundreds  of  spike  returns,  the  probability  of  detecting  spikes 
in  any  individual  waveform  remains  low. 

Equation  (5)  was  used  to  estimate  the  backscattering  cross  section.  For  this  calculation,  the 
volume  backscatter  coefficient  is  assumed  to  be  10"6  m'1  sr'1  [9],  In  Fig.  1,  the  larger  spike 
represents  one  of  the  largest  spikes  detected  during  these  measurements.  In  this  case,  the  ratio 
of  the  power  returns  within  the  bracket  is  18.2  and  the  backscattering  cross  section  is  0.08 
mm2  sr'1.  The  smaller  spike  in  Fig.  1  has  a  ratio  of  0.6  with  a  backscattering  cross  section  of 
0.003  mnrsr"1. 

The  total  density  of  spike  measurements  (with  backscattering  cross  sections  magnitudes 
greater  than  0.003  mm2  sr'1)  was  calculated  and  compared  reported  densities  for  radar 
observations.  Each  lidar  measurement  range  bin  (1.5  m  long)  is  considered  to  be  an 
independent  measurement;  thus  there  are  five  hundred  measurements  for  each  lidar  profile. 
Over  the  measurement  period  interval,  over  ten  thousand  profiles  were  recorded.  Each 
measurement  volume  is  approximately  0.04  m3  and,  with  about  five  million  observations,  the 
total  measurement  volume  is  about  2.0  x  105  m3.  (The  measurements  at  individual  range  bins 
are  assumed  to  independent  of  each  other  which  may  be  true  for  weaker  spikes,  but  not 
always  true  for  stronger  spikes  that  can  impact  more  than  one  range  bin.)  The  target  density 
is  about  5.7  x  10'3  spikes  per  cubic  meter  which  is  significantly  higher  than  the  density  (5  x 
10‘6  m'3)  reported  by  Martin  and  Shapiro  [11]  for  measurements  made  by  a  9.3  GHz  radar. 
This  radar  density  is  the  same  as  the  density  observed  by  the  lidar  for  particulates  with 
backscattering  cross  sections  larger  than  approximately  0.01  mm2  sr'1. 

Figure  3  shows  the  probability  distribution  of  spikes  as  a  function  of  backscattering  cross 
section.  The  probability  distribution  has  a  tail  that  extends  beyond  the  spikes  caused  by 
Gaussian  noise.  The  Gaussian  noise  is  caused  by  an  analog  background  signal  (primarily  sky 
and  electronic  noise)  and  the  tail  that  extends  beyond  the  Gaussian  noise  is  backscatter  from 
particulates.  A  power  law  distribution  was  fit  to  this  tail: 

Pi=a[AP'i)m  ’  (7> 

where  pi  denotes  the  probability  of  finding  a  spike  in  the  ith  measurement  bin,  a  is  the 
amplitude,  A  P'i  is  the  size  of  the  ith  bin,  and  m  is  the  magnitude  of  power  law.  An  uneven, 
geometric-progress  bin  size  is  used  to  generate  the  histogram: 

AP'.=AP'  X‘  ,  (8) 

z  o’  11 

where  APe  is  the  size  of  the  smallest  bin  size  and  X  is  the  step  size  of  the  geometric 
progression.  The  bin  sizes  vary  between  0.003  and  0.08  mm2  sr'1.  The  geometric  bin  size  is 
used  to  improve  the  bin  statistics  for  the  largest  spike  returns.  However  there  are  empty  bins 

for  many  of  the  large  spike  return  values.  We  used  a  least  squares  algorithm  to  fit  the  data. 

The  minimum  bin  size  used  in  the  fit  has  a  spike  magnitude  of  0.003  mm2  sr'1.  The  largest  bin 
size  is  smaller  than  the  first  empty  bin  (i.e.  the  smallest  bin  which  did  not  have  any  recorded 
spikes);  the  maximum  cross  sectional  magnitude  used  is  about  0.02  mm2  sr'1.  The  bins  with 
the  largest  spikes  (which  are  over  0.02  mm2  sr'1)  are  larger  than  the  smallest  empty  bin  and  are 
not  used  in  the  fitting  process. 
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Backscattering  Cross  Section  (mm^r'1) 


Fig.  3.  A  histogram  of  probability  distribution  of  spikes  is  shown  as  a  function  of  the 
backscattering  cross  section.  The  blue  curve  shows  a  Gaussian  distribution  fit  to  cross 
sections  smaller  than  0.003  mm2  sr'1  and  red  curve  shows  a  power  law  fit  to  sizes  larger  than 
0.003  mm2sr'1and  less  than  0.02  mm2sr'1.  The  dotted  line  extensions  to  red  curve  are  provided 
for  visualization  purposes  only  and  not  derived  from  any  fit  beyond  the  above  mentioned  data 
range. 


5.  SUMMARY 

This  paper  reports  on  lidar  measurement  of  spike  returns  made  over  the  Chesapeake  Bay. 
Using  Eq.  (5),  the  range  of  backscattering  cross  sections  is  between  0.003  and  0.08  mm2  sr"1. 
While  these  cross  section  can  not  be  directly  compared  with  the  cross  section  sizes  reported 
by  Martin  and  Shapiro  [11]  for  radar  systems,  the  lidar-observed  measurement  density  (  5.7  x 
10‘3  spikes  per  cubic  meter)  is  larger  than  the  density  reported  by  radar  systems.  The  larger 
measurement  density  can  be  explained  in  part  by  the  sensitivity  of  the  lidar  to  "targets"  with 
smaller  cross  sections. 

The  measurements  in  this  paper  are  derived  solely  from  remote  sensing  measurements.  We 
do  not  have  any  in  situ  or  physical  measurements  of  the  large  particulates,  which  cause  the 
spike  returns,  nor  do  we  know  the  shape  or  composition  of  the  particles.  The  errors  associated 
estimating  the  cross  sections  associated  with  the  spike  returns  are  large  (possibly  an  order  of 
magnitude).  However  the  magnitude  of  spike  returns  is  consistent  with  returns  from 
biological  particulates  (such  as  insects)  or  ultra-giant  aerosols  (UGA)  [12]. 

Published  results  [13]  suggest  that  high  power  millimeter  wave  radars  can  detect  both  hard 
target  and  diffuse  atmospheric  features  associated  with  refractive  structures.  Since  lidar  does 
not  detect  spike  returns  from  diffuse  refractive  structures,  joint  lidar-radar  measurements  of 
spike  returns  provide  a  technique  for  separating  of  radar  spike  returns  from  "spot  angels"  or 
hard  target  scattering  from  more  diffuse  returns  from  atmospheric  refractive  structures.  Lidar 
measurements  can  provide  a  density  measurement  for  larger  aerosols,  which  are  potentially  a 
source  of  atmospheric  clutter  for  millimeter  radars.  Birds  and  insects  can  move  at  speeds  that 
differ  from  the  wind  speed  and,  thus,  represent  a  source  of  systematic  error  for  active- 
microwave  wind  sounders  —  lidar-detected  spike  returns  may  represent  a  way  to  differentiate 
between  moving  targets  (like  birds  and  insects)  and  large  particulates  which  follow  the  air's 
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motion.  If,  indeed,  some  lidar  spike  returns  are  caused  by  UGA,  the  measurements  are 
potentially  significant,  since  UGA  could  significantly  impact  warm  rain  initiation  [14].  Few 
studies  of  UGA  have  been  made  and  lidar  measurements  may  provide  a  technique  for 
tracking  UGA  and  determining  their  atmospheric  residency  times.  For  the  above  reasons,  we 
think  that  lidar  spike  returns  merit  further  study. 

If  future  measurements  are  made,  the  limitations  and  inaccuracies  of  the  observations 
reported  in  this  paper  should  be  reduced.  Calibrated  targets  should  be  used  to  ensure  that  the 
lidar  spike  data  can  be  more  accurately  converted  to  cross  section  size.  In  situ  measurements 
of  the  aerosol  spectra  should  be  made  to  reduce  the  errors  in  estimating  extinction  and 
backscatter.  Joint  lidar  and  in  situ  aerosol  collection  should  be  made  to  provide  detailed 
information  on  the  characteristics  of  the  particulates  which  are  causing  the  spike  returns. 
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